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Abstract 

We propose a dark energy model with a logarithmic cosmological fluid which 
can result in a very small current value of the dark energy density and avoid the 
coincidence problem without much fine-tuning. We construct a couple of dynamical 
models that could realize this dark energy at very low energy in terms of four scalar 
fields quintessence and discuss the current acceleration of the Universe. Numerical 
values can be made to be consistent with the accelerating Universe with adjustment 
of the two parameters of the theory. The potential can be given only in terms of the 
scale factor, but the explicit form at very low energy can be obtained in terms of 
the scalar field to yield of the form V((f>) = exp(— 2^>)(^0 + B). Some discussions 
and the physical implications of this approach are given. 



1 Introduction 



One of the most intriguing discovery of modern cosmology is the acceleration of the Uni- 
verse [H [2] and it is widely believed that dark energy of repulsive nature is causing the 
current acceleration. Many candidates of the dark energy have been proposed [31 EJ E], 
among which the cosmological constant is the most accepted one. Along with yet another 
unidentified constituent of the Universe called dark matter, they compose standard cos- 
mological model, A — CDM |6]. It is remarkable that the observable Universe can be well 
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addressed with the A — CDM model. Still, extreme fine-tuning of the cosmological con- 
stant [7] has been unsatisfactory feature of the model. That is, the current cosmological 
constant must have an the unnaturally small value compared to Planck scale. 

An alternative proposal to explain the dark energy is the quintessence model [HI E] 
in which a scalar field is added as an indispensable component of the universe. In this 
approach, the smallness of the cosmological constant is achieved by a dynamical decay 
of the scalar field energy density. It has the very attractive features of tracking behavior 
and attractor solutions [I] such that galaxy formation is not affected too much by the 
quintessence field and the dark energy becomes dominant only at the late stage of the 
Universe causing the current acceleration. However, in order to achieve the late time dark 
energy dominance, thus providing a possible solution of the coincidence problem [3], the 
theory has to be fine-tuned to a certain extent such that the energy density today must 
be very close to the critical density [TO] . 

In this paper, we propose a dark energy model which can alleviate the fine-tuning 
problem substantially. Suppose a decaying cosmological term [11] according to l/a(t) 2 , 
where a(t) is the scale factor of the Universe [121 EG]- We regard it as the size of the 
Universe. This term decreases with the expansion of the Universe, and the current value is 
A ~ 10~ 122 M 2 , with ao ~ 10 42 Gev _1 , where we assumed that its value at the Planck scale 
is of the order M 2 . Note that energy density of the cosmological constant from this value is 
very close too the critical density p cr ,o — 1CT 122 M 4 , and it has the potential of explaining 
the coincidence problem without fine-tuning. Also, some theoretical background was 
given for such a decay [121 [13]. On the other hand, the conservation of the Einstein 
tensor prevents the cosmological term to be varying in pure gravity. If matter contents are 
included, the varying cosmological constant term disrupts the matter continuity equations, 
which changes the predictions of the standard cosmology in the matter-dominated epoch 
[I2l IT3] . However, if continuity equation is being enforced, 1/a 2 term behaves exactly like 
the curvature constant term and it alone cannot yield an accelerating universe 

It turns out that by adding another cosmological term which varies according to In a/ a 2 
to the original 1/a 2 , accelerating Universe can be realized. To see this in detail, let us first 
check the energy dominance of each epoch. We start from p cr ~ 1CT 122 M P 4 and the matter 
energy density p m = nj^js Pcr, 1 — 0.27 [6]. Let us assume that the energy density of the 
cosmological fluid composed of those two terms was approximately of the order of ~ M p 4 
when the inflation started 0. At the end of the inflation when the scale factor becomes 10 3 
cm with the number of e-foldings given by N, the energy density would have decreased in 

lr rhe cosmological terms considered here decay very fast and cannot be responsible for the inflation 
itself. The inflaton should come from some other source. 
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magnitude by a factor of ~ Ne 2N and becomes of the order of ~ Ne 27V M 4 . We choose 
N ~ 81 which is bigger than the minimum number of e- foldings [10] so that the energy 



density becomes of the order of 
energy density 



10 M p and therefore we propose the following dark 



Pd 



c* m(a/a inf ) 



x 10 _72 M 



p ' 



(a/a inf ) 2 (a/a ml ) 2 _ 

where c* and d* are constants § and a; n f ~ 10 3 cm. Note that the number 10 -72 is a 
dynamical consequence of the inflation and it turns out that no extreme fine-tuning of 
c* and d* are necessary to describe the current accelerating Universe. Comparing with 
Pd,o = 8pcr, 8 — 0.73 [6], we have a relation 



c* (25 In 10) + d* ~ 5 



;i.2) 



Just after the inflation ended at the energy scale ~ 10 Gev, the energy density is 
of the order of 10" 24 M 4 » p D , and p r ~ 10~ 24 /(a/a r>i ) 4 M 4 with a r>i ~ a ini ~ 10 3 cm. 
The normal expansion takes over and the Universe expands by a factor of 10 21 until the 
radiation-matter equality around a ~ 10 cm. When this occurs, p r j ~ p m ^ ~ 
and the matter dominance takes over since the radiation energy density decays faster than 
the matter energy density. In the meantime, the magnitude of po keeps on decreasing 
according to (11.11) . and becomes of the order of ~ 10~ 114 M P 4 . Therefore, the dark energy 
is completely subdominant during this period. 

Then, the matter-dominated epoch began around a ~ 10 _3 ao, ao ~ 10 28 cm. Since, the 
matter energy density decays faster than dark energy, there exist a scale where p m ~ po 
given by 



a 



c* In 



a 



+ 5 



7- 



;i.3) 



So far, there is only one restriction on the numerical values of the parameters c* 
and d* of (11.21) . and a wide range of their values are allowed to fit into the current 
observation. One can impose one more condition by demanding that the Universe has 
begun its acceleration very recently. The acceleration equation with our po and is 
given by 



a 
a 



1 



(a/a ) 2 (a/a ) 3 



Pc 



2 We also assume that the pressure is given by pd = - 
Pd satisfy the continuity equation, po = —3H(pu +pd)- 



c=* ln(a/ ai n f) 
(a/a inf ) a 



c,+d, 
(a/a inf ) a 



x 10" 



2 M p 4 . p D and 
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In the next sections, we will show that this equation can be realized within a couple 
of quintessence models at very low energy. From the above equation, we see that the 
acceleration began around a acc ~ 7/0*00- For example, if we choose c* = 0.54, a acc ~ 
l/2ao and a cq ~ 0.60ao- This would determine ~ —56. Therefore, the transition to 
dark energy dominance occurred very recently. With these values, the dark energy stays 

negative during most of the time until a h ~ 0.25ao and becomes positive at late stage 

of the matter-dominated era. After passing this point a maximum, a max is reached, and 
eventually, it begin to be dominant around a cq . It seems that a priori there is no reason 
for it to stay positive always as long as the total energy density p = po + p r or p = p D + p m 
remains positive. In addition, absolute value of the energy density is very small when it 
stays negative compared to the radiation or matter energy density. Therefore, it should 
not disturb the radiation- or matter-dominated evolution to alter the course of it §. Note 
that the increasing behavior of dark energy, although limited until a max in our case, also 
appears in the phantom model [14] . It is interesting to check that a_+, a max , a acc , a cq all 
happen very recently without too much fine-tuning of the parameters c* and d*. 

In Sec. 2, we construct a quintessence model with four scalar fields which can produce 
the dark energy behavior of (11.11) at very low energy. An explicit form of the potential in 
terms of the scalar field in the scalar-dominated region is given. In sec. 3, we consider a 
generalized quintessence model where the explicit construction can be extended to matter- 
dominated epoch. In Sec. 4, a critical analysis of the generalized quintessence model is 
performed. Sec. 5 contains conclusion and discussions. 

2 Quintessence with four scalar fields 

Let us consider an action of the form (8irG = 1) 



where a a is an 3 component scalar field and a is a parameter which we assume to be 
constant. S m ia a matter action with p m = 0. The spacetime metric tensor is given by 



3 It can be checked that a 1- ~ e s / c *clq independently of the detailed numerical values for p cr and 




with 



X = g^S ab d^ a d v a b , 



(2.5) 




(2.6) 
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where a(t) is the scale factor. With ansatz for the scalar field a a of the form a a = x a [T5] , 
the evolution and the continuity equation for matter are given as follows; 

SH 2 = h 2 + ^ + V + Pm , (2.7) 

2 GT 

-2H = 2 + ^? + p m , (2.8) 



= p m + 3Hp m , (2.9) 

where H = a /a is a Hubble parameter, p m is the matter energy density. The scalar field 
satisfies 

+ 3#0 + V\4>) = , 

where the dot and prime denote partial differentiations with respect to t and 0, respec- 
tively. 

Taking linear combinations of (12.71) and (12. 8p . one gets the acceleration equation 

! = !( v -^H*- (2 ' io > 

In the quintessence model, the kinetic term and matter density decay very fast in the 
above equation (I2.10p , and sole potential dominance at late time is reached which causes 
the acceleration. We open the possibility that the late time acceleration comes from the 
first two terms in the above equation, and anticipating (I2.10p to describe the current 
acceleration at very low energy, we require the following relation; 

v - P = -Tr*** , (2.ii) 

{a/a y 

where A is a parameter which we assume to be the positive. With this, the expansion 
changes from the deceleration for the small value of the scale factor corresponding to the 
matter-dominated epoch to an acceleration for the large value corresponding to the scalar- 
dominated epoch including the kinetic energy density. We will omit a and p CT in what 
follows unless confusion arises. Using the relation (12. lip , the evolution equations (12. 7p 
and (12.81) can be rewritten as follows: 

n 3 A 3a , 
3H = 2 V -M + ^ + »-> < 2 ' 12 > 

-2H = V-4 + ^+A»- (2.13) 

Differentiating equation (12.121) with respect to time and comparing with (12.131) . we obtain 
the following first-order differential equation for the potential 

AA 

3aV (a) + 6V (a) - — = , (2.14) 
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which yields 

4,4 In a B . 

which is precisely of the form (II. ip . From equation (I2.12p . we define (with a = —A/3, 
see eq. (12.181) and below) 

2Alna 3(B - A) 

Comparing with Eq. (Il.ip . we obtain A = c*/2, 3B = 3c*/2 + 26. 

Note that it is difficult to express the scalar field in a closed form from the equation 

U23U, 



+ — — . ( 2 - 17 ) 



'4Alna B — A 

3a? a? 

so the potential 1/ cannot be expressed in terms of the explicitly. Also, positivity of the 
square root in the above expression restricts the applicable range of a which turns out to 

be a > a h which is the same as the positivity of the of (I2.16p . Basically, this comes 

from the imposition of the acceleration condition (12. lip , and this condition restricts the 

dynamically permitted region to a > a |_. A closed form of the potential is viable, if we 

neglect the matter density in the evolution equations, that is, in the scalar-dominated 
epoch. In this case, the ratio of if 2 and 2 is given by 



'9 V-4 



3m fy-^ + r 



(2.18) 



and we can adjust the parameter a = —\A. Then, we have 

0(a) = v^lna + C, (2.19) 
with an integration constant D. This gives the the following expression for the potential 

V{4) = e -W^ViA . (2.20) 

where the term involving the constant C are absorbed into A and B. 

3 Generalized quintessence model 

In this section, we consider a generalized quintessence model given by 

S 2 = J d 4 xy^[^ - l -uj{<P)d^d^ - V(<f>) - al] + S m , (3.21) 
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where we introduced a scalar function lu(4>) which will be arranged for our purpose. The 
last term in the gravity sector is the triplet of scalar fields as before, and we choose the 
same ansatz of the form a a = x a , which solves the equation of motion. The potential 
V{4>) takes the following form 

V{<f>) = e~ 2 * f^ + B) . (3.22) 



^ 3 

It is essentially of the form given by (I2.20p with a redefinition of the scalar field. We 
assume the matter to be a cold dark matter with w m = as before and satisfies the 
continuity equation. The evolution equations are given by 

3H 2 = l ^ 2 + ^ + V + p m , (3.23) 

2 GT 

-2H = u<p 2 + ^ + p m , (3.24) 

or 

along with the continuity equation (I2.9P and the scalar field satisfies 

u<j) + + 3H0 + V = . (3.25) 

Taking linear combinations of (I3.23P and (I3.24p . one gets the acceleration equation 

- = \(v-ujtf)-p m . (3.26) 
a 3 V / 

We choose the same acceleration condition as before 

V - ujtf = 4 , (3.27) 
cr 

with a positive A. Note that Eqs. (13.231) and (13.241) imply (I3.25P and recall p m = j^j^sPcr- 
Therefore, we have two independent dynamical equations and one constraint (I3.27P for 
the three unknown functions to be determined, a(t) or H, u)((f>), and 4>{t). It turns out 
that we can solve the equations with the ansatz 

0(a) = In a, (3.28) 

which reproduces the same forms of the scale factor dependent potential (I2.15P and energy 
density (I2.16p . Inserting this ansatz into (I3.27p . we obtain the following u in terms of <fi 

= 4A0 + 3(5-A) + 2 7 e-^ (3 ' 29) 
Note that for the choice of a = —A/3, lo{4>) > for positive value of the energy density 
defined with the first three terms of eq. (I3.23P as before. One can show that for a < —A/3, 
uj{4>) is always greater than zero for a > a_ + . The above analysis shows that it is possible 
to construct an explicit form of the potential beyond the scalar-dominated region if we 
introduce a generalized quintessence model with an adjustable kinetic function ou(4>). 
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4 Critical Analysis 



In this section, we perform the critical analysis jl] of the evolution equations of the 
previous section and present some numerical result. For convenience, we choose a = —A. 
For other choices, the qualitative feature of the stability does not change as long as 
a < —A/3. Let us introduce the following dimensionless quantities 

X "o^'^3^' (4 - 30) 

with V = V + 3a/ a 2 . Then the equations (I3.23l) - (l3.25p can be written in the following 
form 

= -3x 2 + -x--y, (4.31) 
dN 3 3 y v ' 

dy 1 4 . 



where N = In a and together with a constraint equation 

{& = 1 -*-V' ( 4 - 33 ) 

The critical points of the above system are easily obtained by setting the right-hand sides 
of the above equations (14.311) and (14.321) to zero. The only physically meaningful critical 
points (x c , y c ) of the system are 

(A) : (x e ,y e ) = (0,0), 

(B) : (x c ,y c ) = Q, | 

Point (A) corresponds to the matter-dominated point, whereas point (B) is a scalar- 
dominates one. 

To gain some insight into the property of the critical points, we write the variables 
near the critical points (x c , y c ) in the form 

x = x c + 5x (4.34) 
y = y c + 5y (4.35) 

where 5x and 5y are perturbations around the critical points. From (14.3 lj) and (I4.32p . we 
obtain the linearized equations 

dSx ( 5 \ . 1 . / , x 

w - {- 3 -^ c )5x--Sy, (4.36) 

(i-3y c )8x+(l-3x c )8y, (4.37) 



dN \3 * c \3 
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which can be written by using a matrix M: 





i 



" { ] j - J (438) 

One can study the stability of critical points against perturbations by evaluating the 
eigenvalues of the matrix M. For class (A) corresponding to the matter dominated epoch, 
one has Ai = 1 and A2 = 1, which means that it is an unstable point. For class (B) 
corresponding to the scalar field dominated epoch, we obtain Ai = —1 and A2 = 0. The 
appearance of zero eigenvalue means that the linear perturbation, which leads to the 
matrix (I4.38p . is not adequate and we must consider the higher order perturbations to 
determine whether the considered critical point is stable or not. 
First note that we have 

2x - y > , (4.39) 

for a = —A. Therefore, it is useful to change from the (x,y) to new variables (X,Y) 
defined as 

X = x + 2y , Y = 2x - y , (4.40) 

where the X, Y are positively defined, and then (14.31ft and (14.321) can be rewritten as 

dX "i 1 

- = -Y--X { 3X + e y-B), (4.41) 

% = -\y(3X + 6Y-5), (4.42) 
and corresponding critical points in terms of the new variables are : 

(A) : (A c ,F c ) = (0,0), 

(B) : {X B ,Y C )= 0,0 

In the higher order perturbations, we write the variables near the critical points 
(X c , Y c ) in the form 

X = X c + 5X (1) + 5X (2) + 5X (3) + • • • , 

Y = Y c + 5Y (l) + 5Y (2) + 5Y {3) + ■ ■ • , (4.43) 

where SX^ and SY^ are n-th order perturbations of the variables near the critical 
points. For class (B), the perturbative equations of each order are given by 

— 5X^ = -5X^--5Y^ 
dN 3 
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dN 



d 

dN 



5X^ 



— 8Y< 
dN 



d 

In 



8X^ 



dN 



-hx™ + 25YW) - 8X® - -5Y^ 

5 v ' 3 

-\$Y {X) (8X^ + 28Y^) 

5 5 3 

_^ 5 y(D 5 y(2) _ 3 ( 5X (i) 5 y(2) + OT^Y^) 
5 5 v ' 



Note that the right hand side of the second equation in the first order perturbation 
equation is zero which reflects the fact that one of two eigenvalues is zero and we must 
focus on the next-order equations. 

The solutions of the above linear differential equations are given by 



X(N) = - + 8A 1 + 5B 1 , 
Y{N) = + 8A 2 + 8B 2 , 



(4.44) 
(4.45) 



with 

8A 1 
5A 2 
8B 1 



5Bo 



-V + \SY 2 N- V 3 iV 2 

5Y - 5Y 2 N + 8Y*N 2 + C(e 4 
3 _„„ (\ 



Saoe 



-N 



1 + -e-^ao 



0{e 4 



N SYn 



-2N 



25 



8 at 



25 



-e- N (l 



: 5Y 5aoe 



-N 



e~ N 8a - 2N5Y ) + 0(e 



(4.46) 
(4.47) 
ION) 8a 8Y 

(4.48) 

(4.49) 



where 8Y and 8a$ are the initial values at N = that satisfy the 8Y = 8Y^(0), 
5ao = 8X^(0) + -^8Y^(0), and e is the infinitesimal order parameter for the perturbation. 
If the higher order terms are included and when N becomes very large, 8A\ )2 can be 
expressed in a closed form with 



8A 1 -»• 
8A 2 



8Yn 



3(1 + 8Y N) 
8Y 



8Y N 



SB^O, 
8B 2 -»■ . 



(4.50) 
(4.51) 
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Since the variable Y is a positively defined value, perturbation around the zero point 
must be also positive. In this case, 8A\ and SA2 smoothly go to zero when N goes 
infinity. Therefore, critical point (B) is stable. Numerical result is given in Fig. 1 which 
demonstrates the stability of the scalar-dominated critical point. 




0.0 0.5 1.0 , r 1.5 2.0 2.5 



Figure 1: The flow diagram of the system in terms of X, Y for different initial values. We have 
drawn the physically allowed region with 5Y > 0. The left red dot is the matter-dominated 
point (A), whereas the right red dot is the scalar-dominated point (B). 

5 Conclusion 

We showed that the coincidence problem can be avoided with logarithmic cosmological 
fluid of the form (11. ip . and this can be realized dynamically as a couple of quintessence 
models at very low energy. Among the four scalar fields scalar fields, one of the scalar fields 
play a major role as in the standard quintessence model causing the current acceleration, 
whereas the triplet of scalar fields is not essential in dynamical evolution but it can provide 
the necessary energy density such that the potential is completely integrable as an exact 
expression in the dark-energy dominated era. In the generalized quintessence model, the 
construction was extended to matter-dominated epoch, and critical analysis indicates that 
the scalar-domination is dynamically stable. An analytic expression of the potential in 
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terms of the scalar field is unavailable in each case, but effective field theories could be 
be considered separately at each stage of the evolution of the Universe. We also have 
checked that in such a scheme, the current value of extremely small dark energy density 
can be obtained without much fine-tuning; The constants c* and are only of the order 
10 _1 and 10 2 , and the small value is attributed to its decaying according to essentially as 
1/a 2 . 

Comparing with the quintessence model, the constraint (12. lip brings a crucial differ- 
ence as far as energy dominance is concerned. When the scalar field begins to roll down 
the potential, the initial potential energy is converted into kinetic energy which dominates 
the energy of the scalar field. But the kinetic energy decreases rapidly and the potential 
energy dominance takes over around 1/2 < z < 1, which is responsible for the dark energy 
domination. In our case, at the very low energy eq. (12. lip maintains throughout the evo- 
lution, and it implies that kinetic energy does not decay fast, but it remains comparable 
with potential energy ~ V/2. This would result in a relatively small absolute value for 
the equation of state but detailed comparison with the observational data needs to 
be done. We comment that this constraint is a phenomenological input to conform to the 
current acceleration, and a more theoretical basis is required. 

Another comment is the feature that the dark energy density remains negative during 

most of the time until a h . Even though this behavior does not destroy the accelerating 

Universe, there seems to be no special meaning to a ~ 0.25ao in the evolution. It needs to 
be addressed how the matter-dominated evolution is effected by the small negative energy, 
if any. One could get rid of the negative energy simply by modifying the dark energy 

density into p£>=-RHS of (II. ip for a > a h but it seems rather ad-hoc. If this scale could 

somehow be raised to electro-weak breaking scale, this could endow more flexibility to the 
theory. Adding (In a) 2 /a 2 term to the dark energy density ( ll.ljl is another possible avenue 
to deal with the negative energy. It is likely that this will also extend the applicable range 
of the quintessence model similarly constructed, and might also improve the equation of 
state previously mentioned. 

We conclude with an intriguing property of (11.11) . One can check that that the status 
current accelerating Universe is rather insensitive to the numerical values of a m { ~ 10 3 cm 
and p CT ~ 10~ 122 Mp chosen. If these values are chosen differently, these changes can be 
always reabsorbed into c* and d*. Suppose we change p CT into a~ 1 p cr with some constant 
a. Then, eqs. (11.11) and (11.4j) retain the same form with c* — > ac*, and d* — > ad*. 
Likewise, if a- m i — > /3 _1 ai n f, then, c* — > (3 2 c*, and d* — > /3 2 (c*ln j3 + d*). Since a and (3 
can be at most of the order 1, c* and d* can change in maximum by an order of 2, and 
fine-tuning problem does not arise. And these changes do not alter the energy dominances 
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of radiation epoch and matter epoch. This especially means that a acc can be made in fact 
independent on a in f without much fine-tuning, suggesting that the current status of the 
Universe is not effected too much by the early universe. Similarly, the scaling of the scale 
factor itself can also be absorbed into new definition of c* and d*, which seems to suggest 
a kind of scaling behavior of the dark energy proposed. It would be interesting if some 
theoretical foundation could be given for it. 
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